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Abstract 



We observe a series of discontinuities in the above gap current/voltage 
characteristics of highly transmissive Nb/Si superconductor /semiconductor 
interfaces. The bias at which each switch occurs decreases in quantized steps 
as a function of temperature and magnetic field. The observed step-size is 
the same for both dependencies. We propose that the steps correspond to 
supercurrent quantization through granular Josephson junctions formed near 
the interface. 
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The conductance through a quantum point contact in a two-dimensional electron gas 
(2DEG) is quantized to 2e^/h times the number of transmitted one- dimensional channels 
Analogous behavior is predicted for a small area Josephson junction formed by a narrow 
constriction between two superconductors |@]. In this case, quantization is in the critical 
current, and is given by eAo/h times the number of transmitted one-dimensional channels, 
where Aq is the pairing potential of the superconductor. Recently, Takayanagi et al. 
fabricated the first working superconducting quantum point contacts, using Nb contacts to 
a submicron width InAs 2DEG, modulated by a split-gate. They observed steps in both the 
critical current and the normal conductance as a function of gate bias. However, substantial 
fluctuations were observed in the step height, and the critical current was not quantized at 
the theoretically predicted values. MuUer et al. [|] measured the conductance and critical 
current through a break junction of Nb. They observed steps in the critical current consistent 
with the theoretically predicted value, suggesting that individual conductance channels are 
gradually shut off as the filament is broken. More recent experiments imply, however, 
that such critical current steps might be caused by abrupt atomic rearrangements, which 
change the normal resistance 0. 

This paper presents intriguing new evidence for the quantization of the critical current 
in small area Josephson junctions. We perform detailed measurements of the interface 
resistance of Nb/Si junctions. Below a transition temperature of 1.5 K, we observe a series of 
switches in the above gap DC current /voltage characteristics. As the temperature increases 
towards the transition, or as a magnetic field is applied, the switching bias decreases in a 
series of quantized steps. The observed step-size is equal in both temperature and magnetic 
field dependencies. We propose that the switches occur as we exceed the critical supercurrent 
of Josephson junctions formed between NbsSi particles and the Nb contact. By varying 
the temperature or magnetic field we are able to reduce the junction area, and observe a 
corresponding quantized change in the critical current. 

Fig. 1(a) shows a schematic drawing of our Nb/Si junctions. A Nb contact is made to 
a 53 nm thick single crystal n"'"^ Si layer which has been wafer bonded to 100 nm of Si02 
on a p"*" Si substrate. 2.0 K Hall measurements of the Si layer indicate that the carrier 
density n = 5.3 x lO^^m"^ and the elastic scattering length /g = 10.8 nm. We create a 
highly transmissive interface by first sputter-etching the Si surface at 200 W for 10 minutes 
in an Ar pressure of 3 x 10~^ mbar. This removes a 25 /xm x 25 /im square of the (100) 
Si top-layer, and exposes clean (111) Si side-surfaces. Without breaking the vacuum we 
then deposit 180 nm of Nb, resulting in a total Nb/Si contact area of approximately Ac = 6 
(/zm)^. We measure the critical temperature of our Nb to be = 9.0 K, indicating a high 
quality film. Four-terminal measurements of the Nb/Si junction resistance are made via Ag 
contacts to the Si, and Au contacts to the Nb. The effects we describe below have been 
observed in four different devices, each fabricated using the same method. 

In the inset to Fig. 2, the differential resistance (dV/dl) of a typical junction is plotted 
as a function of applied bias. Measurements at three temperatures (1.6, 1.2 and 0.35 K) are 
shown, normalized with respect to the resistance at high bias. In each trace, there is a broad 
dip in dV/dl corresponding to the Nb gap energy (Aq = ±1.6 meV). Further measurements 
(not shown) demonstrate that the dip survives up to the Nb transition temperature. This 
dV/dl feature has been reported by a number of groups [^|[ and can be attributed to 
Andreev refiection at the Nb/Si interface P,pi][. In the 0.35 K trace we also observe a strong 
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sub-gap resistance dip. (This first appears at temperatures below 1.1 K.) A similar feature 
was reported by Heslinga et al. and attributed to the appearance of a 'proximity induced 
gap,' in the Si layer. More simply, the RSJ (resistively shunted junction) model predicts that 
a sub-gap dip can be observed if a Josephson junction forms in parallel with the interface 
resistance [pH]] . 

In our case, a possible explanation for the sub-gap resistance dip is suggested by the 
temperature dependence of dV/dl at zero bias (not normalized), plotted in the main part 
of Fig. 2. The differential resistance decreases sharply at a temperature Ti = 1.5 K, from 
890 to 480 fl. This drop, which strongly resembles a superconducting transition, occurs at 
the reported transition temperature of NbsSi ||12|. If we assume that the ratio between the 
transition temperature Ti and the gap energy Am-iSi equals that for Nb, then Am^si = 0.25 
meV. This is similar to the observed resistance dip width of ±0.28 mV at r= 0.35 K. These 
results imply that the sub-gap feature is due to Andreev reflection from NbsSi particles that 
form near the interface during fabrication. NbaSi has a high eutectic temperature (645° 
C), but it has been shown that Nb diffuses nonuniformly into Si at our maximum process 
temperature of 300°C (immediately following the sputter-etch), and forms small intermixed 



regions ||T3[. We also measured junctions made without the sputter-etch, and these showed 
no sub-gap dip or resistance transition. 

Figure 3(a) gives results of a four-terminal DC measurement of the Nb/Si interface (see 
inset) at 1.6 and 0.35 K, corresponding to temperatures above and below the transition 
temperature, Ti. At T = 1.6 K, the I-V curve has a weak Schottky barrier dependence. At 
T = 0.35 K, the current rises above the 1.6 K value, as expected from the results of Fig. 
2. As the applied bias increases, the measured current and voltage switch discontinuously 
due to a sudden increase in the interface resistance. (A switch is observed in both the 
current and the voltage across the interface Vm because only the total applied bias Va is 
fixed at each measurement point.) A total of six switches are observed with increasing bias, 
until the above and below transition I-V characteristics are almost identical. In Fig. 3(a), 
we also plot the 0.35 K characteristic measured with a transverse magnetic field of 5 = 
1.0 T (dashed line). The I-V is almost identical to the above transition trace measured at 
T — only a small excess current is observed at low bias. We note that these results are 
completely reproducible, and that non-reproducible noise due to heating is not observed 
until the current > 200 /iA. The maximum power through the device is ~1 /iW, while our 
pumped '^He cryosat is able to maintain a 100 heat load at 0.35 K for 3.5 hours. 

In Fig. 3(b), we plot a series of 26 I-V characteristics for temperatures between 1.6 and 
0.35 K. The lowest lying trace corresponds to r= 1.6 K. The traces above this correspond to 
the temperatures 1.55 K, 1.50 K,. . . , 0.35 K and are offset by 0.25 /iA, 0.50 /lA,. . . , 6.25 /lA, 
respectively. Plotted in this way, the resistance switches reveal a remarkable temperature 
dependence. As the temperature increases from 0.35 K, the bias at which each switch occurs 
decreases discontinuously, in a series of steps. For any particular switch, the step width is 
approximately constant. If we consider the three most clearly defined switches (marked 5*1, 
5*2, and S3), the average width of the first four steps is dVi = 1.60 mV, dV2 = 1.82 mV and 
(iVs = 2.04 mV, with a maximum variation (per switch) of 3 %. The steps can be seen more 
clearly if the switching bias Vn is plotted as a function of temperature. This is shown in the 
inset of Fig. 3(b) where Vn is plotted for the three switches Si, S2 and S3 as a function of the 
reduced temperature T/Ti. The biases are normalized by the average step widths, dVn. In 
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each case, three or more equally spaced steps are observed. As the temperature approaches 
Ti, the switches merge together. 

Further measurements provide more convincing evidence for the quantization of the 
switching bias. Figure 4(a) shows the normalized switching bias as a function of temperature 
taken from observations of the lowest bias switch in a second device. These points are 
compiled from a set of 94 I-V measurements (similar to those shown in Fig. 3(b)) taken at 
equally spaced temperatures between 1.30 and 0.37 K. A total of 12 well-defined steps are 
observed, each containing at least three bias points. The average step width, dV, is 0.42 
mV. We point out that 4 other switches were observed in this device, and each showed a 
similar temperature dependence, although dV is different for each switch. 

We also studied the influence of magnetic field on the switching bias. Figure 4(b) shows 
the normalized switching bias for the switch shown in Fig. 4(a) as a function of magnetic 
field perpendicular to the current path. These points were compiled from a set of 130 I-V 
measurements taken at magnetic fields equally spaced between and 0.33 Tesla (T = 0.37 
K). The data clearly show that the switching bias is also quantized as a function of magnetic 
field. We observe 12 well-defined steps each containing 6 or more i?-field points. The bias 
is again normalized using dV = 0.42 mV. For B > 0.25 T, the steps dissapear, and the 
switching bias decreases contiuously with magnetic field. Most remarkably, a comparison 
between Figs. 4(a) and (b) shows that the quantized step-size is equal to the same value in 
the 5-dependence as in the T-dependence. The background Tand 5-dependences, however, 
are clearly different. Similar results were observed in all 5 switches of this device. 

To understand these results, we now consider the Nb/Si interface in more detail (see 
Fig. 1(b)). Current flows across the interface by electron tunneling through the Si deple- 
tion region. Because of the high doping density, the average inter-dopant spacing (n~^/^ ~ 
2.7 nm) is comparable to the depletion width (3.7 nm) []T1|. In this case, it has been 
shown that the depletion width varies substantially across the contact and significant trans- 
mission occurs only at the most transparent sections of the interface |T^. We calculate 



the approximate area of the transparent regions as follows. The average transmission co- 
efficient per quantum channel is given by tav = Rs/Rj where Rs is the Sharvin resis- 
tance and Rj is the measured junction resistance. From Fig. 2, Rj at 1.6 K = 890 Q, 
and Rs for our contact is 20 mQ so that tav = 2.3 x 10^^. On the other hand, a fit 
of the normalized dV/dl dip at 1.6 K to the BTK formalism gives a transmission coeffi- 
cient for the transparent regions of tA = 0.84 [0. The maximum area of the transparent 
regions is then approximately At = x (tav/tA) = 164 nm^, where Ac is the total con- 
tact area. Below the transition, at T = 1.2 K, the resistance drops to Rj = 480 Q , 
but the size of the normalized dip in dV/dl remains roughly the same (see inset. Fig. 2). 
This means that stays constant, while the transmissive area increases to approximately 
164 nm2 x (890 / 480 fi) = 304 nm^. 

Additional conduction pathways below Ti are created if NbsSi particles form near the 
interface during sputtering. The pathlength for tunneling via the particles is reduced at 
the NbaSi transition temperature, and the effective transmissive area increases. The cross- 
sectional area for the NbsSi particles is approximately equal to the increase in the transmis- 
sive area or f« 140 nm^. As shown in Fig. 1(b), if the normal coherence length C,n reaches 
the NbsSi particles a supercurrent can flow through weak-link Josephson junctions between 
the NbaSi particles and the Nb (Ji and J2 in Fig. 1(b)). 
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Consider now the switches observed in the I-V charactersitics of Fig. 3. The switches 
progressively transform the T < Ti resistance (0.35 K) into the T > Ti resistance (1.6 K) 
with increasing bias. Each switch is thus a partial breakdown of the conduction pathways 
added at Ti. This could be possibly explained by assuming that Joule heating raises the 



temperature of the conduction pathways above Ti. In this case |T6|, the switching bias 
should follow a (Ti — T)^/^ dependence. Our observed switching bias, however, increases 
more slowly at low T than this. In addition, it is unclear how a heating model can explain 
the quantization we observe in Fig. 4. 

Resistance steps have also been reported by Hebard and Vandenberg in the I-V character- 
istics of high-resistance granular lead films |[T^ . These were later explained as the switching 
of rows of underdamped Josephson junctions formed between grains, lying perpendicular to 
the current path |T^. A switch occurs as a row moves from a zero- voltage state, to one 
in which a vortex undergoes a phase-slip process. The switches we observe might possibly 
be due to a similar process ocurring among granular Josephson junction arrays within the 
depletion region. Typically, though, the separation between switches in a voltage biased 
Josephson-j unction array are independent of small magnetic field, and as the temperature 
is increased, the switches simply disappear [0,|1^] — quantization is not observed. 

We might assume, then, that a group of Josephson junction pathways are associated 
with each particular switch. As the coherence length in the depletion region decreases with 
increasing T or B, individual conduction channels are cut off, resulting in discrete steps. 
But since the observed step spacing is quantized for a particular switch, this would require 
an identical contribution from each random pathway. 

We instead propose that the switches are due to the critical current transition of individ- 
ual nanometer scale Josephson junctions formed between NbaSi particles and the Nb contact. 
Consider a simple circuit model for the Nb/Si contact (Fig. 1(c)). Transport through the 
nth silicide particle is given by a Josephson junction J„ in series with a NbsSi/Si interface 
resistance In the superconducting state, the Josephson junction resistance is zero, and 
the current through J„ is simply /„ = V/Rn- The switching bias is thus proportional to the 
current carried by the Josephson junction at the critical current transition. 

If the area of the Josephson junction is small enough, the switching bias should be quan- 
tized to V = Ne Aq/ Rnh, where N is the number of one-dimensional transmitted channels, 
and eAo/ Rnh is equal to the quantized step size dV. In our case, N = kpA/An, where A is 
the cross-sectional area of the junction. To produce the observed quantization, we assume 
that A decreases as a function of B and T. Since the normal coherence length ^„ decreases 
with increasing B and T, this is reasonable. The number of modes A^in each junction (Fig. 3 
and 4) is simply equal to the normalized switching bias V/ dV. From the maximum number 
of modes, we can calculate the approximate cross-sectional area An of the particles corre- 
sponding to the three switches 5*1, S2, and 5*3 in Fig. 3. This gives Ai = 56 nm^, A2 = 84 
nm^ and ^3 = 131 nm^, roughly corresponding to the cross-sectional area determined from 
our simple estimate above (140 nm^), and giving an average particle radius on the order of 

We expect that, in the limit of large N, the T and B dependences in Fig. 4 should be 
described by Ic(T, B) for a large area point contact. The envelope of our quantized T- 
dependence appears similar to IdT) for a large area device [|20|, while the envelope of the 
quantized 5-dependence resembles the low field behavior of a screened Josephson junction 
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||TT| . The field required to insert a single flux quantum in such small junctions (~ 100 
nm^) is larger than the critical fleld of Nb, so that no oscillatory dependence is expected. A 
detailed understanding would require a rigorous analysis of Ic{T, B) in the crossover between 
macroscopic and nano-scale junctions. 

In conclusion, we observe evidence that the maximum supercurrent through Josephson 
junction pathways formed near a Nb/Si interface is quantized, in agreement with theoret- 
ical predictions. This could provide a new and straightforward way to study supercurrent 
transport in nanometer-scale Josephson junctions. 

The authors thank F. Hekking, T. Klapwijk, M. de Jong, and D. Williams for valuable 
discussions and S. Newcomb for providing TEM analysis. 
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FIGURES 



FIG. 1. (a) Schematic drawing of our Nb/Si junctions, (b) A representative region of the 
Nb/Si interface. NbaSi particles, which form during Nb deposition, provide additional tunneling 
pathways. At low enough temperatures the normal coherence length ^„ extends to create Josephson 
junctions Ji and J2. (c) Circuit model of the Nb/Si interface. The NbsSi Josephson junctions J„ 
have normal state resistances R'^ and lie in series with the Nb3Si/Si interface resistances The 
in-parallel Schottky barrier resistance is R. 

FIG. 2. Zero-bias differential resistance of a highly transmissive Nb/Si junction verses temper- 
ature. Ti is the NbsSi transition temperature. Inset: The normalized differential resistance verses 
bias at 1.6, 1.2 and 0.35 K. The 1.2 and 1.6 K traces are offset by 0.1 and 0.2, respectively. The 
gap widths for Nb and NbsSi are indicated. 

FIG. 3. (a) Four-terminal measurement of the DC current and voltage across the Nb/Si in- 
terface at T = 0.35 K and T = 1.6 K. Each trace shows measurements made at 1000 equally 
spaced values of the applied bias Va (see inset). The dashed line shows the characteristics for B = 
1.0 T and T = 0.35 K. (b) DC current and voltage characteristics plotted for 26 equally spaced 
temperatures between 0.35 K and 1.6 K. The traces are offset for clarity. Inset The normalized 
switching bias as a function of reduced temperature for the three switches indicated. 

FIG. 4. The normalized switching bias as a function of (a) reduced temperature and (b) 
magnetic field for a typical resistance switch. 
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This figure "Fig 1. GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9707071v2 



This figure "Fig2.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9707071v2 



This figure "Fig3.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9707071v2 



This figure "Fig4.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9707071v2 



